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Abstract 

Results of the study of the e + e~ — > 7r°7 process with SND detector at 
VEPP-2M collider in the c.m.s. energy range y/s = 0.60 — 0.97 GeV arc 
presented. Using 36513 selected events corresponding to a total integrated 
luminosity of 3.4p6 _1 the e + e~ — * 7r°7 cross section was measured. The 
energy dependence of the cross section was analyzed in the framework of 
the vector meson dominance model. The data are well described by a sum 
of 4>,w,p — > 7r°7 decay contributions with measured decay probabilities: 
Br{cu -> 7T°7) = (9.34±0.15±0.31) % and Br{p° -> n°j) = (5.15±1.16± 
0.73) x 10~ 4 . The p — uj relative interference phase is ip p ^ = — 10. 2 ±6.5 ± 
2.5°. 

1 Introduction 

Cross section of the e + e~ — > 7r°7 process at the c.m.s. energies y/s = 0.60 — 
0.97 GeV within the framework of the vector meson dominance model is de- 
termined by radiative decays of light vector mesons p°(770), w(782), </>(1020). 
These decays belong to the class of magnetic dipole transitions and represent 
major interest for study of quark structure of vector mesons and for tests of 
low-energy models of strong interactions, such as non-relativistic quark model, 
effective potential models, etc. PHI El E] EJ- Study of this process allows to 
improve accuracy of the parameters of the p°,uj —> ir°~f decays. 

The only previous measurement of the decay p° — * 7r°7 was carried out 
by ND detector 0: Br(p° -> 7r° 7 ) = (7.9 ± 2.0) x 10~ 4 . This result agrees 
with the PDG value for isotopically complementary channel: i3r(p ± — * 7r ± 7) = 
(4.5 ± 0.5) x 10 -4 [Hil- The u> — > 7r°7 decay was studied in several experiments 



H3 DHH21 USUI El- The world average decay probability Br{u -► 7r 7) = 

(8.7 ±0.4)% ma. 

In this work we present the study of the e + e~ — > 7r°7 process with SND 
detector at VEPP-2M collider. 

2 Detector and experiment 

The SND detector consists of an electromagnetic calorimeter, tracking 
and muon systems. The main part of the detector is a three-layer spherical 
electromagnetic calorimeter consisting of 1600 Nal(Tl) crystals. Total thick- 
ness of the calorimeter for the particles flying from the interaction point is 
13.4 Xq. Total solid angle is 90% • 4-7T. The energy resolution of the calorime- 
ter for photons is cr E /E w 4.2%/i?(GeV) 1 / 4 , the angular resolution is cr^g » 
0.82°/ y/E(GeV) 0.63°. 

The experiment was carried out at VEPP-2M collider |T7] with SND detector 
|16j . The data were collected in March - July, 1998 J2| a t 30 points in the 
energy range y/s = 0.60 — 0.97 GeV. The total integrated luminosity used for 
the analysis was ~ 3.4 pb _1 . The beam energy determination was based on 
measured magnetic field in the bending magnets and beam revolution frequency 
in the collider. The error of the center of mass energy determination consists 
of two parts: 0.1 MeV — relative accuracy of energy setting for each energy 
point and 0.2 MeV — general energy scale bias common for all points within 
the experiment. 

3 Data analysis 

In this work the process e + e~ — ► tt°j was studied in the three-photon final 
state. The main sources of background are QED processes e + e~ — ► 37 and 
e+e~ — > 27 with extra photons of the machine background. Other possible 
sources are process e + e~ — > 777, and cosmic background. 

3.1 Events selection 

For an event to be recorded the SND first level trigger (FLT) required at least 
two clusters of hit crystals in the calorimeter and no signals in neither tracking 
nor muon systems. The FLT threshold on calorimeter energy deposition changed 
with the beam energy, but was always below 0Ay/s. 

The reconstructed events were first put through primary selection which 
required at least 3 neutral and no charged particles, total energy deposition 
E to t > 0.65-^3, total momentum measured by the calorimeter P tot < 0.3-y/s, 
polar angles of the two highest energy photons 36° < #1,2 < 144°, the polar 
angle of the third photon (descending order in energy) 27° < 63 < 153°, and 
the energy deposition of this photon E p a > 0.ly/s. These conditions select 
three-photon events, suppress machine background and two photon annihilation 
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events with additional background clusters in the calorimeter. As a result 52415 
events were selected for further analysis. 

In order to improve energy and angular resolution for photons the selected 
events were kinematically fitted with total energy and momentum conservation 
constraints. The fit results are the value of x§ 7 of the hypothesis and fitted 
kinematical parameters of the photons. The kinematic fit improves resolution 
in invariant mass of photon pairs from ir° decays from 11.2 MeV to 8.6 MeV 
(Fig©. 
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Figure 1: Invariant mass distribution for photons pairs in e + e~ — > 7r°7 events 
before (left) and after (right) kinematic fitting. Solid line — MC simulation, 
points — data (y/s = 782 MeV). 

The x§ 7 distribution is shown in Fig0 For additional suppression of cosmic 
and machine backgrounds we required x§ 7 < 20. This cut also implicitly limits 
maximum energy of initial state radiation (ISR) photons in the process under 
study (Fig|2Jl. In order to suppress 2-photon annihilation background, each 
selected event was kinematically fitted to e + e _ — + 27 hypothesis and restriction 
on calculated x% was applied: x% — X% < 0- The x§ 7 — xl 7 distribution is 
shown in Fig [2] 

The only significant backgrounds to the process under study remaining after 
described above cuts are e + e _ — > 777 and QED 37 annihilation. In the latter 
process all kinematically allowed combinations of the photon energies and angles 
are present, so this background cannot be completely eliminated and must be 
subtracted. To this end all events, which passed primary selection and kinematic 
fit cuts, were divided into two classes: events with 108 MeV < m 77 < 162 MeV 
were assigned to class A, the rest — to class B. Here m 77 is an invariant mass 
of a photon pair after kinematic fitting (Fig©. The total number of the selected 
class A events is 36513. The fraction of e + e~ — > 777 background in this class 
is less than 0.1%. For the class B this fraction is up to 5%. For calculation 
of integrated luminosity the special sample of e + e~ — > 27 events (class C) 
was selected using following criteria: no charged particles, at least two neutral 
particles, energy depositions for two most energetic photons £^1,2 > 0.3y/~s, 
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Figure 2: Distributions over x§ 7 (left) and x| 7 — X27 (right) for class A selection. 
Solid line — MC simulation, points — experiment (y/s = 782 MeV). 

their polar angles 36° < #1,2 < 144°, azimuth acollinearity Aipi 2 < 10°, polar 
acollinearity A#i 2 < 25°, event does not belong to classes A or B. It is necessary 
to note significant contribution of e + e~ — > 7r°7 events to the class C (up to 10% 
at w(782) resonance). 



3.2 Cross section parameterization 

The e + e~ — * 7r°7 cross section in the framework of VDM can be parameterized 
as follows EH EHJ: 
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Here g 7 y and gvn°-y are coupling constants, my is the V resonance mass, ry(s) 
is the energy-dependent width of the resonance taking into account processes 
with branching ratios larger than 1%, A nonres represents possible non-resonant 
contribution. Using following formulas for coupling constants: 
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, Eq.QJ can be transformed to the form more suitable for data 
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approximation: 



q(s) 3 



cr 7r o 7 (s) = — g-^- \A p07 ,o 7 (s) + A W7r o 7 (s) + A V7r7 (s) + a 7r o 7 | , (6) 



m v r y /v(a) / m|- 
£>y(s) y g(mf,) d 

where a 7r o 7 is a non-resonant contribution. We used two different models for de- 
scription of interference phase between p, ui — > 7r°7 decay amplitudes. For model 
with energy- independent interference phases, f Ptlp = e lVp , / w = 1. In this case 
the (fp is expected to be 0° for pure p and u> isotopic states. Electromagnetic 
p-to mixing leads to nonzero value of ip p . It can be estimated from B(u> — ► 27r): 
(f p as —13°. The second model is based on mixed propagator approach |27l I28| : 



r„(») = !+«(»)• ™ + , (9) 
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<•„(») = i-e(.)-[f^r ■ r^rl- 

where H puJ is p — uj mixing self-energy. 

Detection efficiency for e + e~ — > 7r°7 process depends not only on c.m.s. 
energy -y/s but also on energy of extra photons emitted by initial particles E r . 
The detection efficiency £ r (y / s, E r ) was determined by Monte Carlo simulation 
with ISR taken into account. The energy and angular distributions of ISR 
photons were generated according to Refs. |21U22| . The dependence of detection 
efficiency on E r approximated by a smooth function is shown in Fig|3]for ^fs = 
782 MeV. The noticeable peak in detection efficiency near ISR kinematic limit 
corresponds to the case when 7r°7 invariant mass is close to m n o and ISR photon 
is emitted at the large angle and detected. The effective threshold on the ISR 
photon energy SE r (s) is determined by x 2 restriction and can be defined as a 
width at half maximum of the e r ( v / s, E r ). At y/s = 782 MeV SE r w 64.4. The 
SE r dependence on \fs is shown in Fig [21 

The visible cross section of the e + e~ — > 7r°7 process was calculated as |2(J) : 



(T V i s (s) = / £ r (Vs, —t- )F(x, s)er((l - x)s)dx (12) 

Jo z 

where <j(s) is cross section, the function F(x, s) is electron "radiator" function 
|21| . For data presentation we used traditional form: 

(Jvis(s) = e(\/i) • P(Vs) ■ cr(s) 
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Figure 3: Detection efficiency as function of the extra photons energy e r (s/s, E r ) 
for y/s = 782 MeV (left) and effective extra photons energy threshold SE r (right) 
as function of y/s. 



where £(\/s) and (3(y/s) are defined as: 

e(y/7s) = e r (y/7s,0) (13) 

0{V - s) = Jo ^ e r (yfs, x -^)F(x,s)a((l~x) S )dx 

s r {y/s, 0) • a(s) 

For simulation of the background process e + e~ — > 37 (QED) the lowest- 
order formulas from were used. Visible cross section calculated using Monte 
Carlo simulation was corrected for higher order loop diagrams and soft photons 
emission [21] using 5E r as an upper limit of soft photons energy. The correc- 
tion varied in the range of 0.915 — 0.925. We expect that the accuracy of the 
calculated e + e~ — > 37 visible cross section is not worse than 2%. 

For simulation of the process e + e~ — > 27 (QED) used for luminosity deter- 
mination the formula from taking into account additional photon emission 
was used. The accuracy of the visible cross section determination is estimated 
as 1%. 

3.3 Data approximation 

The FIT package [20] was used for data fitting. The fitting was done by means of 
maximum likelihood method on all three data sets (classes A,B,and C) simul- 
taneously. Expected number of events in the i-th energy point was calculated 



as: 
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Visible hadronic cross sections were calculated according to Ea. H12|) and cor- 
rected for beam energy spread. Because the e + e~ — ► ir°~f process gives notice- 
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able contribution to events of luminosity process e + e~ — > 2j, the integrated 
luminosity (ILi) was recalculated on every iteration step of the minimization. 
The accuracy of determination of the c.m.s. energy is worse than the accuracy 
of the w-meson mass value. Therefore we introduced the possible energy scale 
bias AE as a free parameter. Other fit parameters were er W7r o 7 , cr p7r o 7 , a 7r o 7 , ip p 
or U pi0 depending on approach to phase factor (Eqs.Q-ljHJl) calculation, and 
k^y. The k^j parameter is a ratio of measured and calculated cross sections. 
Parameters of the process e + e~ — * <fr — > ir°j were taken from |29| : 

cr 0w o 7 = 5.12 ± 0.39 nb (15) 
W = 158 ±11° (16) 

For other cross section parameters the world average values were used. 
Data were approximated in four models: 

1. cr p o 7r o 7 = 0, a 7r o 7 = 0, 

2. CTpO„.o 7 and tp p o are free parameters, a w o 7 = 0, 

3. cr p o T o 7 is a free parameter, a w o 7 = 0, n paJ is a free parameter, 

4. CT p o 7r o 7 is a free parameter, a 7r o 7 is a free real parameter, IL puJ calculated 
from B(u> — > 2ir). 

Obtained energy scale bias (AE) is —0.34 ± 0.08 MeV for all models. It is 
consistent with our expectations. Found value of is 98.7±1.3% with 2^ = || 
shows good agreement between calculated and measured cross sections of QED 
37 annihilation. For background subtraction we used the measured cross section. 
Other obtained parameters are listed in the Table □ 



Table 1: The fitted cross section parameters for different models. Only statisti- 
cal errors are shown. 





<W 7 ' n k 


cr p o OT o 7 , nb 




n pw ,MeV 2 


Re a 7r o 7 ,nbi 


X 2 /N 


1 


176.6±1.4 












81/28 


2 


155.8±2.7 


0.58±0.13 


-10.2±6.5 







21.6/26 


4 


155.9±2.7 


0.56±0.13 


-9.9±6.5 1) 


-2819±1841 





21.9/26 


3 


156.8±2.8 


0.51±0.13 


-12.8±1.1 1) 


-3676±303 2) 


-0.13±0.13 


20.7/25 



^ calculated using Eqs.(JEjHdJ 

2 ) derived from B(u — > 2it) 



Large \ 2 value for the first model shows that e + e~ — > 7r°7 cross section 
cannot be described only by u and </> decays contribution. The second model 
corresponds to an energy independent p — lo interference phase. Obtained value 
of this phase (—10.2 ± 6.5)° is in agreement with expectation from electro- 
magnetic p — lo mixing tp p = (—12.8 ± 1.1)°. Therefore the last two fits were 
performed in mixed propagator approach (Eq.JSJ)). Mixing self-energy H puJ was 
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taken as a free parameter for the model 3 and calculated from the world average 
B(u! — ► 27r) for the model 4. The model 4 was used to estimate the contribution 
from higher vector resonances p', u>', which was introduced as a pure real param- 
eter a v Oj. The fitted value of a 7r o 7 is compatible with zero. All the models 2 — 4 
describe the experimental data equally well. 

3.4 Systematic errors. 

Systematic error contributions for obtained cross section parameters are sum- 
marized in the Tabled 

Systematic error of luminosity determination originates mostly from inaccu- 
racy in cross section calculation (1%) and uncertainty of detection efficiency of 
the luminosity process, which was estimated using different angle and acollinear- 
ity selection cuts. The total error of the integrated luminosity determination 
was ~ 2 — 3%. 

Primary selection efficiency depends on simple kinematic cuts and is inde- 
pendent of c.m.s. energy. Thus, its systematic error emerging from simulation 
inaccuracy was studied by comparison of simulated and experimental event dis- 
tributions at cj-resonance peak, where backgrounds are negligible. Systematic 
error of primary selection efficiency does not exceed 1.5%. 

Table 2: Contributions to the systematic errors of the cross section parameters 



Source 


CW 7 


CT p 7r 7 


<P 


Integrated luminosity 


2.0% 


3.1% 


8% 


Three photons selection efficiency 


1.5% 


1.5% 


4% 


Final selection efficiency 


1.6% 


5.6% 


22% 


Additional clusters 


0.3% 


1% 


2% 


PDG table errors 


0.1% 


4% 


4% 


Total (no model error) 


3.0% 


9.0% 


24% 



The machine background changed with c.m.s. energy. To study its influ- 
ence on detection efficiency we merged recorded experimental background events 
with the simulated ones. Comparison of detection efficiencies obtained by sim- 
ulation with and without merging machine background gives an estimate of the 
detection efficiency error from this source, not exceeding 0.5%. 

The final class A selection conditions contain cuts in invariant masses and 
complex kinematic parameters x| 7 j xi 7 - Dependences of their efficiencies on 
c.m.s. energy and ISR photon energy in experiment and simulation may differ. 
In order to evaluate systematic error coming from this source, approximations 
were done with different cuts in these parameters. 

Substantial systematic error contributions to the a p a T .a 1 and ip come from 
inaccuracy of PDG data, mostly from the r w uncertainty. 
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Table 3: e + e~ — > 7r°7 cross section. is a c.m.s. energy spread, 7X is an 
integrated luminosity, N is a number of events, Nb g is an estimated number of 
background events, e 7r o 7 is detection efficiency (Ea. l|13|l ) of the process e + e~ — ► 
7r° 7 , /3 7r o 7 is a factor taking into account radiative correction (Eq.(|2|) and 
energy spread. Energy (y/s) is corrected according to fitted AE, its error is 
0.08 MeV. The first error of the cross section <j v o~ is statistical, the second one 



is systematic. 



y/s, McV 


SE, McV 




JV 










599.52 


0.14 


39.90 ± 0.30 


60 


46.0 


0.315 


0.908 


1.23 ±0.77 ±0.19 


629.51 


0.15 


46.09 ± 0.33 


62 


44.8 


0.316 


0.904 


1.31 ±0.68 ±0.43 


659.52 


0.16 


40.02 ± 0.33 


37 


32.0 


0.313 


0.899 


0.45 ±0.63 ±0.17 


689.56 


0.19 


48.31 ± 0.38 


48 


33.0 


0.316 


0.893 


1.10 ±0.59 ±0.30 


719.51 


0.18 


58.43 ± 0.43 


69 


36.0 


0.323 


0.886 


1.97 ±0.56 ±0.18 


749.50 


0.20 


50.90 ±0.42 


84 


26.8 


0.317 


0.866 


4.09 ±0.73 ±0.28 


759.50 


0.20 


41.88 ± 0.39 


107 


20.7 


0.316 


0.846 


7.71 ±1.02 ±0.42 


763.50 


0.21 


38.80 ±0.38 


124 


18.7 


0.317 


0.834 


10.27 ±1.19 ±0.47 


769.50 


0.21 


43.60 ±0.40 


234 


20.3 


0.319 


0.812 


18.95 ± 1.45 ± 1.08 


773.50 


0.21 


62.77 ±0.48 


531 


28.7 


0.319 


0.794 


31.62 ± 1.51 ± 1.45 


777.50 


0.21 


76.73 ± 0.53 


1544 


35.0 


0.319 


0.776 


79.60 ±2.13 ±2.15 


778.71 


0.22 


6.88 ±0.16 


162 


3.2 


0.319 


0.772 


93.89 ±8.13 ±3.82 


779.48 


0.24 


43.39 ± 0.39 


1282 


20.0 


0.319 


0.770 


118.49 ±3.46 ±2.94 


780.59 


0.23 


132.36 ±0.68 


4989 


61.4 


0.319 


0.772 


151.36 ±2.20 ±4.40 


781.63 


0.24 


351.62 ± 1.10 


15259 


164.1 


0.319 


0.779 


172.82 ±1.43 ±5.15 


782.52 


0.21 


81.11 ± 0.53 


3523 


37.8 


0.319 


0.793 


169.83 ±2.94 ±4.76 


783.51 


0.21 


74.90 ± 0.51 


3150 


34.7 


0.319 


0.816 


159.53 ±2.93 ±5.46 


785.51 


0.22 


73.73 ± 0.51 


2391 


33.3 


0.320 


0.883 


113.15 ±2.39 ±3.18 


789.50 


0.22 


56.91 ± 0.46 


930 


24.6 


0.320 


1.044 


47.61 ± 1.66 ± 1.43 


793.49 


0.23 


53.03 ± 0.45 


456 


22.4 


0.320 


1.201 


21.26 ±1.10 ±0.80 


799.49 


0.23 


51.86 ± 0.45 


285 


21.3 


0.320 


1.411 


11.27 ±0.77 ±0.31 


809.49 


0.25 


65.73 ± 0.52 


189 


25.8 


0.318 


1.660 


4.70 ±0.43 ±0.13 


819.49 


0.24 


115.74 ±0.70 


233 


43.7 


0.318 


1.775 


2.90 ±0.25 ±0.19 


839.47 


0.25 


144.83 ± 0.80 


179 


52.6 


0.320 


1.711 


1.59 ±0.18 ±0.07 


879.45 


0.27 


170.26 ±0.91 


100 


50.4 


0.318 


1.269 


0.72 ±0.16 ±0.17 


919.43 


0.32 


327.70 ± 1.32 


137 


77.7 


0.319 


1.048 


0.54 ±0.12 ±0.04 


939.45 


0.30 


291.22 ± 1.28 


99 


63.0 


0.318 


1.007 


0.39 ±0.12 ±0.09 


949.45 


0.29 


259.10 ± 1.22 


86 


53.9 


0.317 


0.993 


0.39 ±0.13 ±0.06 


957.45 


0.29 


241.63 ± 1.18 


79 


48.8 


0.317 


0.984 


0.40 ±0.13 ±0.06 


969.46 


0.30 


245.65 ± 1.21 


84 


47.6 


0.319 


0.969 


0.48 ±0.13 ±0.10 



4 Results 

Our final results are based on model 2 approximation. Differences in approxi- 
mation results for models 2-4 were considered as model error contributions to 
total systematic errors. As a result we present: 

cr e + e -^_ 7r o 7 = (155.8 ± 2.7 ± 4.8) nb, (17) 
(T e +e-^p°^n°-y = (0.58 ± 0.13 ± 0.08) nb (18) 
<f> P u = (-10.2 ±6.5 ±2.5) degrees (19) 

Detailed point by point listing of the measured e + e _ — > ir°j cross section is 
presented in Table |2I Systematic error of the experimental cross section is 



9 



determined by systematic errors of integrated luminosity, detection efficiency 
radiative correction, and background subtraction. It is worth mentioning that 
systematic errors for different c.m.s. energy points are highly correlated. Mea- 
sured cross section and data from [2H] are also plotted in FigQ] 




500 600 700 800 900 1000 1100 

y/s, MeV 

Figure 4: e + e~ — * 7r°7 cross section. Solid line depicts the cross section in the 
model 2, dashed line — in the model 1. Data from ND experiment are grouped 
by energies and shifted taking into account current world average value of the 
w-meson mass. 

Decay parameters expressed in terms of probabilities and partial widths are: 



Br(w - 


* 7T°7) ■ B(u -> e+e~) 


= (6.50 ± 0.11 ± 0.20) x 10" 6 , 


(20) 


Br(p - 


-> 7T 7 ) ■ B(p -> e+e") 


= (2.34 ± 0.53 ± 0.33) x 10" 8 , 


(21) 




Br(uj -> 7T°7) 


= (9.34 ±0.15 ±0.31)%, 


(22) 




Br{p° -> 7T 7 ) 


= (5.15 ± 1.16 ± 0.73) x 10" 4 , 


(23) 




— >7r°7 


= 788 ± 12±27keV, 


(24) 




r p o^ Tr o 7 


= 77±17±llkeV 


(25) 



Obtained results statistically agree with previous measurements. The partial 
width I\,o _>„.<>_, is in a good agreement with the world average V p ±^^± 1 . Phc- 
nomenological estimates using various models |3 H Ej do not contradict 
our result. 

The ratio of the partial widths of the u>,p — > ir°j decays required by strict 
SU(3) symmetry [2] is equal to 9.47, which agrees with our measurement: 



7 = 10.3 ± 2.5 ± 1.4, (26) 
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5 Conclusions 



The most accurate measurement of e + e — > 7r°7 cross section is performed at 
the c.m.s. energy region of 0.60 - 0.97 GeV at the VEPP-2M collider with the 
SND detector. At present experimental accuracy level this cross section is well 
described by vector meson dominance model taking into account transitions 
4>,uj,P — > 7r°7. In this model the cross sections of the processes e + e~ — > u> — > 
7r°7 and e + e~ — ► p — ► 7r°7 at corresponding meson masses are measured. 
Partial widths, their rations and decay probabilities of corresponding decays 
were evaluated. Results are presented in Eas. l|17[l - i|26[l and in the Table|3J 

Measured values of the u, p° — > 7r°7 decay parameters are consistent with 
earlier experimental results. Partial width of the p° — > 7r°7 decay is in a good 
agreement with that of p — > 7r ± 7 decays. These values also do not contradict 
various phenomenological estimations. Obtained value of the p — to interference 
phase could be well explained by electromagnetic p — lu mixing. Our results 
have higher accuracy than the world average for p° — > 7r°7 and u) — > ir 0r y. 
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